Interactions between calcium-induced arrhythmia triggers and the electrophysiological-anatomical substrate underlying the induction of atrial fibrillation.

Atrial fibrillation (AF) is the most common cardiac arrhythmia and is sustained by spontaneous focal excitations and re-entry. Spontaneous electrical firing in the pulmonary vein (PV) sleeves is implicated in AF generation. The aim of this simulation study was to identify the mechanisms determining the localisation of AF triggers in the PVs and their contribution to the genesis of AF. A novel biophysical model of the canine atria was used that integrates stochastic, spontaneous subcellular Ca2+ release events (SCRE) with regional electrophysiological heterogeneity in ionic properties and a detailed three-dimensional model of atrial anatomy, microarchitecture and patchy fibrosis. Simulations highlighted the importance of the smaller inward rectifier potassium current (IK1 ) in PV cells compared to the surrounding atria, which enabled SCRE more readily to result in delayed-afterdepolarisations that induced triggered activity. There was a leftward shift in the dependence of the probability of triggered activity on sarcoplasmic reticulum Ca2+ load. This feature was accentuated in 3D tissue compared to single cells (Δ half-maximal [Ca2+ ]SR  = 58 µM vs. 22 µM). In 3D atria incorporating electrical heterogeneity, excitations preferentially emerged from the PV region. These triggered focal excitations resulted in transient re-entry in the left atrium. Addition of fibrotic patches promoted localised emergence of focal excitations and wavebreaks that had a more substantial impact on generating AF-like patterns than the PVs. Thus, a reduced IK1 , less negative resting membrane potential, and fibrosis-induced changes of the electrotonic load all contribute to the emergence of complex excitation patterns from spontaneous focal triggers. KEY POINTS: Focal excitations in the atria are most commonly associated with the pulmonary veins, but the mechanisms for this localisation are yet to be elucidated. We applied a multi-scale computational modelling approach to elucidate the mechanisms underlying such localisations. Myocytes in the pulmonary vein region of the atria have a less negative resting membrane potential and reduced time-independent potassium current; we demonstrate that both of these factors promote triggered activity in single cells and tissues. The less negative resting membrane potential also contributes to heterogeneous inactivation of the fast sodium current, which can enable re-entrant-like excitation patterns to emerge without traditional conduction block.

Patchy fibrosis promotes trigger-substrate interactions that both generate and maintain atrial fibrillation.

Fibrosis has been mechanistically linked to arrhythmogenesis in multiple cardiovascular conditions, including atrial fibrillation (AF). Previous studies have demonstrated that fibrosis can create functional barriers to conduction which may promote excitation wavebreak and the generation of re-entry, while also acting to pin re-entrant excitation in stable rotors during AF. However, few studies have investigated the role of fibrosis in the generation of AF triggers in detail. We apply our in-house computational framework to study the impact of fibrosis on the generation of AF triggers and trigger-substrate interactions in two- and three-dimensional atrial tissue models. Our models include a reduced and efficient description of stochastic, spontaneous cellular triggers as well as a simple model of heterogeneous inter-cellular coupling. Our results demonstrate that fibrosis promotes the emergence of focal excitations, primarily through reducing the electrotonic load on individual fibre strands. This enables excitation to robustly initiate within these single strands before spreading to neighbouring strands and inducing a full tissue focal excitation. Enhanced conduction block can allow trigger-substrate interactions that result in the emergence of complex, re-entrant excitation patterns. This study provides new insight into the mechanisms by which fibrosis promotes the triggers and substrate necessary to induce and sustain arrhythmia.

Exploring interpretability in deep learning prediction of successful ablation therapy for atrial fibrillation.

Background: Radiofrequency catheter ablation (RFCA) therapy is the first-line treatment for atrial fibrillation (AF), the most common type of cardiac arrhythmia globally. However, the procedure currently has low success rates in dealing with persistent AF, with a reoccurrence rate of ∼50% post-ablation. Therefore, deep learning (DL) has increasingly been applied to improve RFCA treatment for AF. However, for a clinician to trust the prediction of a DL model, its decision process needs to be interpretable and have biomedical relevance. Aim: This study explores interpretability in DL prediction of successful RFCA therapy for AF and evaluates if pro-arrhythmogenic regions in the left atrium (LA) were used in its decision process. Methods: AF and its termination by RFCA have been simulated in MRI-derived 2D LA tissue models with segmented fibrotic regions (n = 187). Three ablation strategies were applied for each LA model: pulmonary vein isolation (PVI), fibrosis-based ablation (FIBRO) and a rotor-based ablation (ROTOR). The DL model was trained to predict the success of each RFCA strategy for each LA model. Three feature attribution (FA) map methods were then used to investigate interpretability of the DL model: GradCAM, Occlusions and LIME. Results: The developed DL model had an AUC (area under the receiver operating characteristic curve) of 0.78 ± 0.04 for predicting the success of the PVI strategy, 0.92 ± 0.02 for FIBRO and 0.77 ± 0.02 for ROTOR. GradCAM had the highest percentage of informative regions in the FA maps (62% for FIBRO and 71% for ROTOR) that coincided with the successful RFCA lesions known from the 2D LA simulations, but unseen by the DL model. Moreover, GradCAM had the smallest coincidence of informative regions of the FA maps with non-arrhythmogenic regions (25% for FIBRO and 27% for ROTOR). Conclusion: The most informative regions of the FA maps coincided with pro-arrhythmogenic regions, suggesting that the DL model leveraged structural features of MRI images to identify such regions and make its prediction. In the future, this technique could provide a clinician with a trustworthy decision support tool.

Imaging and biophysical modelling of thrombogenic mechanisms in atrial fibrillation and stroke.

Atrial fibrillation (AF) underlies almost one third of all ischaemic strokes, with the left atrial appendage (LAA) identified as the primary thromboembolic source. Current stroke risk stratification approaches, such as the CHA2DS2-VASc score, rely mostly on clinical comorbidities, rather than thrombogenic mechanisms such as blood stasis, hypercoagulability and endothelial dysfunction-known as Virchow's triad. While detection of AF-related thrombi is possible using established cardiac imaging techniques, such as transoesophageal echocardiography, there is a growing need to reliably assess AF-patient thrombogenicity prior to thrombus formation. Over the past decade, cardiac imaging and image-based biophysical modelling have emerged as powerful tools for reproducing the mechanisms of thrombogenesis. Clinical imaging modalities such as cardiac computed tomography, magnetic resonance and echocardiographic techniques can measure blood flow velocities and identify LA fibrosis (an indicator of endothelial dysfunction), but imaging remains limited in its ability to assess blood coagulation dynamics. In-silico cardiac modelling tools-such as computational fluid dynamics for blood flow, reaction-diffusion-convection equations to mimic the coagulation cascade, and surrogate flow metrics associated with endothelial damage-have grown in prevalence and advanced mechanistic understanding of thrombogenesis. However, neither technique alone can fully elucidate thrombogenicity in AF. In future, combining cardiac imaging with in-silico modelling and integrating machine learning approaches for rapid results directly from imaging data will require development under a rigorous framework of verification and clinical validation, but may pave the way towards enhanced personalised stroke risk stratification in the growing population of AF patients. This Review will focus on the significant progress in these fields.

Corrigendum: Reinforcement Learning to Improve Image-Guidance of Ablation Therapy for Atrial Fibrillation.

[This corrects the article DOI: 10.3389/fphys.2021.733139.].

Reinforcement Learning to Improve Image-Guidance of Ablation Therapy for Atrial Fibrillation.

Atrial fibrillation (AF) is the most common cardiac arrhythmia and currently affects more than 650,000 people in the United Kingdom alone. Catheter ablation (CA) is the only AF treatment with a long-term curative effect as it involves destroying arrhythmogenic tissue in the atria. However, its success rate is suboptimal, approximately 50% after a 2-year follow-up, and this high AF recurrence rate warrants significant improvements. Image-guidance of CA procedures have shown clinical promise, enabling the identification of key patient anatomical and pathological (such as fibrosis) features of atrial tissue, which require ablation. However, the latter approach still suffers from a lack of functional information and the need to interpret structures in the images by a clinician. Deep learning plays an increasingly important role in biomedicine, facilitating efficient diagnosis and treatment of clinical problems. This study applies deep reinforcement learning in combination with patient imaging (to provide structural information of the atria) and image-based modelling (to provide functional information) to design patient-specific CA strategies to guide clinicians and improve treatment success rates. To achieve this, patient-specific 2D left atrial (LA) models were derived from late-gadolinium enhancement (LGE) MRI scans of AF patients and were used to simulate patient-specific AF scenarios. Then a reinforcement Q-learning algorithm was created, where an ablating agent moved around the 2D LA, applying CA lesions to terminate AF and learning through feedback imposed by a reward policy. The agent achieved 84% success rate in terminating AF during training and 72% success rate in testing. Finally, AF recurrence rate was measured by attempting to re-initiate AF in the 2D atrial models after CA with 11% recurrence showing a great improvement on the existing therapies. Thus, reinforcement Q-learning algorithms can predict successful CA strategies from patient MRI data and help to improve the patient-specific guidance of CA therapy.

Toward Patient-Specific Prediction of Ablation Strategies for Atrial Fibrillation Using Deep Learning.

Atrial fibrillation (AF) is a common cardiac arrhythmia that affects 1% of the population worldwide and is associated with high levels of morbidity and mortality. Catheter ablation (CA) has become one of the first line treatments for AF, but its success rates are suboptimal, especially in the case of persistent AF. Computational approaches have shown promise in predicting the CA strategy using simulations of atrial models, as well as applying deep learning to atrial images. We propose a novel approach that combines image-based computational modelling of the atria with deep learning classifiers trained on patient-specific atrial models, which can be used to assist in CA therapy selection. Therefore, we trained a deep convolutional neural network (CNN) using a combination of (i) 122 atrial tissue images obtained by unfolding patient LGE-MRI datasets, (ii) 157 additional synthetic images derived from the patient data to enhance the training dataset, and (iii) the outcomes of 558 CA simulations to terminate several AF scenarios in the corresponding image-based atrial models. Four CNN classifiers were trained on this patient-specific dataset balanced using several techniques to predict three common CA strategies from the patient atrial images: pulmonary vein isolation (PVI), rotor-based ablation (Rotor) and fibrosis-based ablation (Fibro). The training accuracy for these classifiers ranged from 96.22 to 97.69%, while the validation accuracy was from 78.68 to 86.50%. After training, the classifiers were applied to predict CA strategies for an unseen holdout test set of atrial images, and the results were compared to outcomes of the respective image-based simulations. The highest success rate was observed in the correct prediction of the Rotor and Fibro strategies (100%), whereas the PVI class was predicted in 33.33% of the cases. In conclusion, this study provides a proof-of-concept that deep neural networks can learn from patient-specific MRI datasets and image-derived models of AF, providing a novel technology to assist in tailoring CA therapy to a patient.

Identifying locations of re-entrant drivers from patient-specific distribution of fibrosis in the left atrium.

Clinical evidence suggests a link between fibrosis in the left atrium (LA) and atrial fibrillation (AF), the most common sustained arrhythmia. Image-derived fibrosis is increasingly used for patient stratification and therapy guidance. However, locations of re-entrant drivers (RDs) sustaining AF are unknown and therapy success rates remain suboptimal. This study used image-derived LA models to explore the dynamics of RD stabilization in fibrotic regions and generate maps of RD locations. LA models with patient-specific geometry and fibrosis distribution were derived from late gadolinium enhanced magnetic resonance imaging of 6 AF patients. In each model, RDs were initiated at multiple locations, and their trajectories were tracked and overlaid on the LA fibrosis distributions to identify the most likely regions where the RDs stabilized. The simulations showed that the RD dynamics were strongly influenced by the amount and spatial distribution of fibrosis. In patients with fibrosis burden greater than 25%, RDs anchored to specific locations near large fibrotic patches. In patients with fibrosis burden below 25%, RDs either moved near small fibrotic patches or anchored to anatomical features. The patient-specific maps of RD locations showed that areas that harboured the RDs were much smaller than the entire fibrotic areas, indicating potential targets for ablation therapy. Ablating the predicted locations and connecting them to the existing pulmonary vein ablation lesions was the most effective in-silico ablation strategy.

Prolonged ursodeoxycholic acid administration reduces acute ischaemia-induced arrhythmias in adult rat hearts.

Acute myocardial ischaemia and reperfusion (I-R) are major causes of ventricular arrhythmias in patients with a history of coronary artery disease. Ursodeoxycholic acid (UDCA) has previously been shown to be antiarrhythmic in fetal hearts. This study was performed to investigate if UDCA protects against ischaemia-induced and reperfusion-induced arrhythmias in the adult myocardium, and compares the effect of acute (perfusion only) versus prolonged (2 weeks pre-treatment plus perfusion) UDCA administration. Langendorff-perfused adult Sprague-Dawley rat hearts were subjected to acute regional ischaemia by ligation of the left anterior descending artery (10 min), followed by reperfusion (2 min), and arrhythmia incidence quantified. Prolonged UDCA administration reduced the incidence of acute ischaemia-induced arrhythmias (p = 0.028), with a reduction in number of ventricular ectopic beats during the ischaemic phase compared with acute treatment (10 ± 3 vs 58 ± 15, p = 0.036). No antiarrhythmic effect was observed in the acute UDCA administration group. Neither acute nor prolonged UDCA treatment altered the incidence of reperfusion arrhythmias. The antiarrhythmic effect of UDCA may be partially mediated by an increase in cardiac wavelength, due to the attenuation of conduction velocity slowing (p = 0.03), and the preservation of Connexin43 phosphorylation during acute ischaemia (p = 0.0027). The potential antiarrhythmic effects of prolonged UDCA administration merit further investigation.

Image-Based Computational Evaluation of the Effects of Atrial Wall Thickness and Fibrosis on Re-entrant Drivers for Atrial Fibrillation.

Introduction: Catheter ablation (CA) is a common treatment for atrial fibrillation (AF), but the knowledge of optimal ablation sites, and hence clinical outcomes, are suboptimal. Increasing evidence suggest that ablation strategies based on patient-specific substrates information, such as distributions of fibrosis and atrial wall thickness (AWT), may be used to improve therapy. We hypothesized that competing influences of large AWT gradients and fibrotic patches on conductive properties of atrial tissue can determine locations of re-entrant drivers (RDs) sustaining AF. Methods: Two sets of models were used: (1) a simple model of 3D atrial tissue slab with a step change in AWT and a synthetic fibrosis patch, and (2) 3D models based on patient-specific right atrial (RA) and left atrial (LA) geometries. The latter were obtained from four healthy volunteers and two AF patients, respectively, using magnetic resonance imaging (MRI). A synthetic fibrotic patch was added in the RA and fibrosis distributions in the LA were obtained from gadolinium-enhanced MRI of the same patients. In all models, 3D geometry was combined with the Fenton-Karma atrial cell model to simulate RDs. Results: In the slab, RDs drifted toward, and then along the AWT step. However, with additional fibrosis, the RDs were localized in regions between the step and fibrosis. In the RA, RDs drifted toward and anchored to a large AWT gradient between the crista terminalis (CT) region and the surrounding atrial wall. Without such a gradient, RDs drifted toward the superior vena cava (SVC) or the tricuspid valve (TSV). With additional fibrosis, RDs initiated away from the CT anchored to the fibrotic patch, whereas RDs initiated close to the CT region remained localized between the two structures. In the LA, AWT was more uniform and RDs drifted toward the pulmonary veins (PVs). However, with additional fibrotic patches, RDs either anchored to them or multiplied. Conclusion: In the RA, RD locations are determined by both fibrosis and AWT gradients at the CT region. In the LA, they are determined by fibrosis due to the absence of large AWT gradients. These results elucidate mechanisms behind the stabilization of RDs sustaining AF and can help guide ablation therapy.

Algorithms for left atrial wall segmentation and thickness - Evaluation on an open-source CT and MRI image database.

Structural changes to the wall of the left atrium are known to occur with conditions that predispose to Atrial fibrillation. Imaging studies have demonstrated that these changes may be detected non-invasively. An important indicator of this structural change is the wall's thickness. Present studies have commonly measured the wall thickness at few discrete locations. Dense measurements with computer algorithms may be possible on cardiac scans of Computed Tomography (CT) and Magnetic Resonance Imaging (MRI). The task is challenging as the atrial wall is a thin tissue and the imaging resolution is a limiting factor. It is unclear how accurate algorithms may get and how they compare in this new emerging area. We approached this problem of comparability with the Segmentation of Left Atrial Wall for Thickness (SLAWT) challenge organised in conjunction with MICCAI 2016 conference. This manuscript presents the algorithms that had participated and evaluation strategies for comparing them on the challenge image database that is now open-source. The image database consisted of cardiac CT (n=10) and MRI (n=10) of healthy and diseased subjects. A total of 6 algorithms were evaluated with different metrics, with 3 algorithms in each modality. Segmentation of the wall with algorithms was found to be feasible in both modalities. There was generally a lack of accuracy in the algorithms and inter-rater differences showed that algorithms could do better. Benchmarks were determined and algorithms were ranked to allow future algorithms to be ranked alongside the state-of-the-art techniques presented in this work. A mean atlas was also constructed from both modalities to illustrate the variation in thickness within this small cohort.

Modeling Left Atrial Flow, Energy, Blood Heating Distribution in Response to Catheter Ablation Therapy.

Introduction: Atrial fibrillation (AF) is a widespread cardiac arrhythmia that commonly affects the left atrium (LA), causing it to quiver instead of contracting effectively. This behavior is triggered by abnormal electrical impulses at a specific site in the atrial wall. Catheter ablation (CA) treatment consists of isolating this driver site by burning the surrounding tissue to restore sinus rhythm (SR). However, evidence suggests that CA can concur to the formation of blood clots by promoting coagulation near the heat source and in regions with low flow velocity and blood stagnation. Methods: A patient-specific modeling workflow was created and applied to simulate thermal-fluid dynamics in two patients pre- and post-CA. Each model was personalized based on pre- and post-CA imaging datasets. The wall motion and anatomy were derived from SSFP Cine MRI data, while the trans-valvular flow was based on Doppler ultrasound data. The temperature distribution in the blood was modeled using a modified Pennes bioheat equation implemented in a finite-element based Navier-Stokes solver. Blood particles were also classified based on their residence time in the LA using a particle-tracking algorithm. Results: SR simulations showed multiple short-lived vortices with an average blood velocity of 0.2-0.22 m/s. In contrast, AF patients presented a slower vortex and stagnant flow in the LA appendage, with the average blood velocity reduced to 0.08-0.14 m/s. Restoration of SR also increased the blood kinetic energy and the viscous dissipation due to the presence of multiple vortices. Particle tracking showed a dramatic decrease in the percentage of blood remaining in the LA for longer than one cycle after CA (65.9 vs. 43.3% in patient A and 62.2 vs. 54.8% in patient B). Maximum temperatures of 76° and 58°C were observed when CA was performed near the appendage and in a pulmonary vein, respectively. Conclusion: This computational study presents novel models to elucidate relations between catheter temperature, patient-specific atrial anatomy and blood velocity, and predict how they change from SR to AF. The models can quantify blood flow in critical regions, including residence times and temperature distribution for different catheter positions, providing a basis for quantifying stroke risks.

Novel MRI Technique Enables Non-Invasive Measurement of Atrial Wall Thickness.

Knowledge of atrial wall thickness (AWT) has the potential to provide important information for patient stratification and the planning of interventions in atrial arrhythmias. To date, information about AWT has only been acquired in post-mortem or poor-contrast computed tomography (CT) studies, providing limited coverage and highly variable estimates of AWT. We present a novel contrast agent-free MRI sequence for imaging AWT and use it to create personalized AWT maps and a biatrial atlas. A novel black-blood phase-sensitive inversion recovery protocol was used to image ten volunteers and, as proof of concept, two atrial fibrillation patients. Both atria were manually segmented to create subject-specific AWT maps using an average of nearest neighbors approach. These were then registered non-linearly to generate an AWT atlas. AWT was 2.4 ± 0.7 and 2.7 ± 0.7 mm in the left and right atria, respectively, in good agreement with post-mortem and CT data, where available. AWT was 2.6 ± 0.7 mm in the left atrium of a patient without structural heart disease, similar to that of volunteers. In a patient with structural heart disease, the AWT was increased to 3.1 ± 1.3 mm. We successfully designed an MRI protocol to non-invasively measure AWT and create the first whole-atria AWT atlas. The atlas can be used as a reference to study alterations in thickness caused by atrial pathology. The protocol can be used to acquire personalized AWT maps in a clinical setting and assist in the treatment of atrial arrhythmias.

Novel Computational Analysis of Left Atrial Anatomy Improves Prediction of Atrial Fibrillation Recurrence after Ablation.

The left atrium (LA) can change in size and shape due to atrial fibrillation (AF)-induced remodeling. These alterations can be linked to poorer outcomes of AF ablation. In this study, we propose a novel comprehensive computational analysis of LA anatomy to identify what features of LA shape can optimally predict post-ablation AF recurrence. To this end, we construct smooth 3D geometrical models from the segmentation of the LA blood pool captured in pre-procedural MR images. We first apply this methodology to characterize the LA anatomy of 144 AF patients and build a statistical shape model that includes the most salient variations in shape across this cohort. We then perform a discriminant analysis to optimally distinguish between recurrent and non-recurrent patients. From this analysis, we propose a new shape metric called vertical asymmetry, which measures the imbalance of size along the anterior to posterior direction between the superior and inferior left atrial hemispheres. Vertical asymmetry was found, in combination with LA sphericity, to be the best predictor of post-ablation recurrence at both 12 and 24 months (area under the ROC curve: 0.71 and 0.68, respectively) outperforming other shape markers and any of their combinations. We also found that model-derived shape metrics, such as the anterior-posterior radius, were better predictors than equivalent metrics taken directly from MRI or echocardiography, suggesting that the proposed approach leads to a reduction of the impact of data artifacts and noise. This novel methodology contributes to an improved characterization of LA organ remodeling and the reported findings have the potential to improve patient selection and risk stratification for catheter ablations in AF.

Atrial Heterogeneity Generates Re-entrant Substrate during Atrial Fibrillation and Anti-arrhythmic Drug Action: Mechanistic Insights from Canine Atrial Models.

Anti-arrhythmic drug therapy is a frontline treatment for atrial fibrillation (AF), but its success rates are highly variable. This is due to incomplete understanding of the mechanisms of action of specific drugs on the atrial substrate at different stages of AF progression. We aimed to elucidate the role of cellular, tissue and organ level atrial heterogeneities in the generation of a re-entrant substrate during AF progression, and their modulation by the acute action of selected anti-arrhythmic drugs. To explore the complex cell-to-organ mechanisms, a detailed biophysical models of the entire 3D canine atria was developed. The model incorporated atrial geometry and fibre orientation from high-resolution micro-computed tomography, region-specific atrial cell electrophysiology and the effects of progressive AF-induced remodelling. The actions of multi-channel class III anti-arrhythmic agents vernakalant and amiodarone were introduced in the model by inhibiting appropriate ionic channel currents according to experimentally reported concentration-response relationships. AF was initiated by applied ectopic pacing in the pulmonary veins, which led to the generation of localized sustained re-entrant waves (rotors), followed by progressive wave breakdown and rotor multiplication in both atria. The simulated AF scenarios were in agreement with observations in canine models and patients. The 3D atrial simulations revealed that a re-entrant substrate was typically provided by tissue regions of high heterogeneity of action potential duration (APD). Amiodarone increased atrial APD and reduced APD heterogeneity and was more effective in terminating AF than vernakalant, which increased both APD and APD dispersion. In summary, the initiation and sustenance of rotors in AF is linked to atrial APD heterogeneity and APD reduction due to progressive remodelling. Our results suggest that anti-arrhythmic strategies that increase atrial APD without increasing its dispersion are effective in terminating AF.

Slow Conduction in the Border Zones of Patchy Fibrosis Stabilizes the Drivers for Atrial Fibrillation: Insights from Multi-Scale Human Atrial Modeling.

Introduction: The genesis of atrial fibrillation (AF) and success of AF ablation therapy have been strongly linked with atrial fibrosis. Increasing evidence suggests that patient-specific distributions of fibrosis may determine the locations of electrical drivers (rotors) sustaining AF, but the underlying mechanisms are incompletely understood. This study aims to elucidate a missing mechanistic link between patient-specific fibrosis distributions and AF drivers. Methods: 3D atrial models integrated human atrial geometry, rule-based fiber orientation, region-specific electrophysiology, and AF-induced ionic remodeling. A novel detailed model for an atrial fibroblast was developed, and effects of myocyte-fibroblast (M-F) coupling were explored at single-cell, 1D tissue and 3D atria levels. Left atrial LGE MRI datasets from 3 chronic AF patients were segmented to provide the patient-specific distributions of fibrosis. The data was non-linearly registered and mapped to the 3D atria model. Six distinctive fibrosis levels (0-healthy tissue, 5-dense fibrosis) were identified based on LGE MRI intensity and modeled as progressively increasing M-F coupling and decreasing atrial tissue coupling. Uniform 3D atrial model with diffuse (level 2) fibrosis was considered for comparison. Results: In single cells and tissue, the largest effect of atrial M-F coupling was on the myocyte resting membrane potential, leading to partial inactivation of sodium current and reduction of conduction velocity (CV). In the 3D atria, further to the M-F coupling, effects of fibrosis on tissue coupling greatly reduce atrial CV. AF was initiated by fast pacing in each 3D model with either uniform or patient-specific fibrosis. High variation in fibrosis distributions between the models resulted in varying complexity of AF, with several drivers emerging. In the diffuse fibrosis models, waves randomly meandered through the atria, whereas in each the patient-specific models, rotors stabilized in fibrotic regions. The rotors propagated slowly around the border zones of patchy fibrosis (levels 3-4), failing to spread into inner areas of dense fibrosis. Conclusion: Rotors stabilize in the border zones of patchy fibrosis in 3D atria, where slow conduction enable the development of circuits within relatively small regions. Our results can provide a mechanistic explanation for the clinical efficacy of ablation around fibrotic regions.

Towards patient-specific modelling of lesion formation during radiofrequency catheter ablation for atrial fibrillation.

Radiofrequency catheter ablation procedures are a first-line method of clinical treatment for atrial fibrillation. However, they suffer from suboptimal success rates and are also prone to potentially serious adverse effects. These limitations can be at least partially attributed to the inter- and intra- patient variations in atrial wall thickness, and could be mitigated by patient-specific approaches to the procedure. In this study, a modelling approach to optimising ablation procedures in subject-specific 3D atrial geometries was applied. The approach enabled the evaluation of optimal ablation times to create lesions for a given wall thickness measured from MRI. A nonliner relationship was revealed between the thickness and catheter contact time required for fully transmural lesions. Hence, our approach based on MRI reconstruction of the atrial wall combined with subject-specific modelling of ablation can provide useful information for improving clinical procedures.

Evolution and pharmacological modulation of the arrhythmogenic wave dynamics in canine pulmonary vein model.

AIMS: Atrial fibrillation (AF), the commonest cardiac arrhythmia, has been strongly linked with arrhythmogenic sources near the pulmonary veins (PVs), but underlying mechanisms are not fully understood. We aim to study the generation and sustenance of wave sources in a model of the PV tissue. METHODS AND RESULTS: A previously developed biophysically detailed three-dimensional canine atrial model is applied. Effects of AF-induced electrical remodelling are introduced based on published experimental data, as changes of ion channel currents (ICaL, IK1, Ito, and IKur), the action potential (AP) and cell-to-cell coupling levels. Pharmacological effects are introduced by blocking specific ion channel currents. A combination of electrical heterogeneity (AP tissue gradients of 5-12 ms) and anisotropy (conduction velocities of 0.75-1.25 and 0.21-0.31 m/s along and transverse to atrial fibres) can results in the generation of wave breaks in the PV region. However, a long wavelength (171 mm) prevents the wave breaks from developing into re-entry. Electrical remodelling leads to decreases in the AP duration, conduction velocity and wavelength (to 49 mm), such that re-entry becomes sustained. Pharmacological effects on the tissue heterogeneity and vulnerability (to wave breaks and re-entry) are quantified to show that drugs that increase the wavelength and stop re-entry (IK1 and IKur blockers) can also increase the heterogeneity (AP gradients of 26-27 ms) and the likelihood of wave breaks. CONCLUSION: Biophysical modelling reveals large conduction block areas near the PVs, which are due to discontinuous fibre arrangement enhanced by electrical heterogeneity. Vulnerability to re-entry in such areas can be modulated by pharmacological interventions.